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ABSTRACT 

^ '. We present the results of the first imaging X-ray observations in the 

t;^^ ■ 0.4-10.0 keV bandpass of the nearby starburst galaxies M82 and NGC 

^ ; 253. The ASCA spectra of both M82 and NGC 253 are complex with 

0^\ strong line emission from O, Ne, Fe, Mg, S, and Si, allowing elemental 

abundances to be estimated in the X-ray band for the first time in these 
^ ■ sources. Two components are required to fit the spectra of both galaxies, 

with a "soft" component well described by a thermal model with a tem- 
'O ! perature of ~ 10^^^ K and a "hard" component well described by either 

a thermal model {Thard ~ 10^ K) or a power-law model (F ~ 1.8 — 2.0). 
We find that different models (with different continua) yield absolute 
abundances that differ by more than an order of magnitude, while rela- 
tive abundances are more robust and suggest an underabundance of Fe 
(inferred from the Fe-L complex) relative to a-burning elements. Dust 
depletion may be responsible for the low relative abundance of Fe in- 
ferred from the soft-component fits. We set the most reliable limits to 
date on the Fe-K emission line equivalent width (EW), with an upper 
limit in NGC 253 of EW < 180 eV while a line at ~ 6.6 keV is marginally 
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significant (at the 3cr level) with EW - 100 eV in M82. The low Fe-K 
line emission EW limits suggest either significantly sub-solar abundances 
in the material producing the hard component (if thermal) or that there 
is a significant non-thermal or non-equilibrium contribution to the hard 
component. Most of the soft flux (which originates mostly within the 
central kpc of M82 and NGC 253) is consistent with starburst models of 
supernovae-heated ISM and, to a lesser extent, starburst-driven super- 
wind emission and the direct emission from supernova (SN). The hard 
component in both galaxies may have some contribution from ~ 10^ 
K superwind emission or individual SN, although most of the emission 
probably originates in point sources (most likely blackhole candidates 
or mini-AGN with MEdd > 3 — 2OM0) and, possibly, inverse-Compton 
scattering of IR photons. The similarity of the spectral X-ray character- 
istics of NGC 253 and M82 to some LINERs and low-luminosity AGN 
suggests a hnk between AGN and starbursts (e.g., both may contain an 
accretion-driven emission component). 

Subject headings: galaxies: starburst - galaxies: abundances - galaxies: 
individual (NGC 253, M82) - X-rays: galaxies 



1. Introduction 

M82 and NGC 253 are two of the near- 
est "starburst" galaxies (we assume distances 
of 3.6 and 2.5 Mpc, respectively; Freed- 
man et al. 1994; Davidge & Prichett 1990) 
where extensive star formation in the < 1 
kpc-scale nuclear regions (Rieke et al. 1980) 
has resulted in a spectacular X-ray morphol- 
ogy, with emission concentrated both in the 
nucleus (at levels of one to two orders of 
magnitudes higher than typically seen in nor- 
mal spirals) and extended along their mi- 
nor axes to scales of up to ~ 10 kpc. M82 
and NGC 253 have been extensively stud- 
ied in X-rays, with data from the Einstein 
HRI (Watson, Stanger, & Griffiths 1984; Fab- 
biano & Trinchieri 1984), ROSATBRl (Breg- 
man, Schulman, & Tomisaka 1995) and EX- 
OSAT LE (Schaaf et al. 1989) revealing the 
X-ray morphology below ~ 4 keV, showing 
that much of the nuclear emission is concen- 
trated in multiple point sources. Data from 
the Einstein IPC (Fabbiano 1988), EXOSAT 
ME (Schaaf et al. 1989), Ginga (Ohashi et 
al. 1990; Tsuru et al. 1990) and BBXRT (Pe- 
tre 1994) suggest a complex spectrum with at 
least two components, although the statistics, 
spectral resolution, and/or bandpass were some- 
what limited with these instruments. These 
studies have suggested that the nuclear X- 
ray emission is due to an extended compo- 
nent probably consisting of the emission of 
unresolved supernovae (SN) and SN-heated 
ISM, and individual point sources which are 
likely to be very luminous X-ray binaries and 
young supernova (SN) or possibly even a hid- 
den AGN. These studies also suggest that 
the "plumes" breaking out of their galac- 
tic disks are evidence for SN-driven "super- 
winds" (Heckman, Armus & Miley 1990) based 
on the observed temperatures and luminosi- 



ties. Inverse-Compton (IC) scattering of the 
large IR flux {Ljr ~ lO^^L©; Telesco 1988) by 
relativistic electrons has also been suggested 
as a significant contributor to the X-ray fiux 
(Rieke et al. 1980; Schaaf et al. 1989; Moran 
& Lehnert 1996), a view (possibly) supported 
by the detection of MeV emission in the vicin- 
ity of NGC 253 (Bhattacharya et al. 1994). 
However, the expected intensity of the IC 
component is subject to large uncertainties 
depending on model assumptions (see §4). 

The significant spatial and spectral com- 
plexity suggests that proper deconvolution of 
the spectra would require good spatial and 
spectral resolution. Here we present high- 
quality spectra of M82 and NGC 253 in the 
0.4-10.0 keV bandpass obtained by the X-ray 
satellite ASCA which carries the first imag- 
ing, medium-resolution X-ray spectrometers. 
Although the spatial resolution of ASCA is 
inferior to that of ROSAT, ASCA can con- 
strain the source of the X-ray emission to 
within ~ 1', while previous detectors with 
significant effective area above 2 keV were es- 
sentially non-imaging (although BBXRT al- 
lowed for some crude spatial analysis). We 
show below that the spectra are indeed com- 
plex, showing strong line emission and requir- 
ing multiple components. 

2. The ASCA Data 

M82 was observed by ASCA (Tanaka, In- 
oue, & Holt 1994) on May 19, 1993 and NGC 
253 was observed on June 10, 1993. Briefly, 
ASCA consists of two solid-state imaging spec- 
trometers (SIS; hereafter SO and SI) with an 
approximate bandpass of 0.4-10.0 keV and 
two gas imaging spectrometers (CIS; here- 
after S2 and S3) with an approximate band- 
pass of 0.8-10.0 keV. The SIS observations 
were done in 4-ccd mode although only the 



data from SO chips 1 and 2 and SI chips 
and 3 (containing most of the source counts) 
were used in the present analysis. Times of 
high background were excluded and hot pix- 
els (in the SIS data) were removed. In M82 
and NGC 253, most of the extended emis- 
sion is concentrated within ~ 1' of the nu- 
cleus (Fabbiano & Trinchieri 1984; Watson, 
Stanger & Griffiths 1984; Bregman, Schul- 
man, & Tomisaka 1995). In NGC 253, several 
point sources are clearly present at distances 
of up to several arcminutes from the nucleus 
(see below). Accordingly, source counts were 
extracted from a 6' radius circle centered on 
the nuclear source in each instrument (source 
regions sizes of 3' and 4' are typically cho- 
sen for the SIS and GIS, respectively, in the 
case of a single point source observation). The 
background was determined from the remain- 
ing counts in the chips beyond 6. '5 for the SIS, 
from a 8-13' annulus for M82 GIS data and 
from an offset 5' circle for the NGC 253 GIS 
data. The observed M82 count rates were 0.6- 
1.0 counts s~^ (per detector) among the four 
detectors with net exposure times of ~ 17- 
20 ks and the background comprising 4-6% 
of the total. The observed NGC 253 count 
rates were 0.1-0.3 counts s"^ among the four 
detectors with net exposure times of ~ 27-30 
ks and the background comprising 12-19% of 
the total. 

ROSAT PSPC images of M82 and NGC 
253 show that their extended halos are soft 
(i.e., kT ^ 0.5 keV), contributing negligibly in 
the ~ 0.5 — 2.0 keV bandpass compared to the 
X-ray flux originating within the 6' source re- 
gions described above. Nevertheless, it is pos- 
sible that the background regions used in ex- 
tracting the ASCA background spectra may 
be contaminated by the scattering of source 
photons or, possibly, some (relatively hard) 



non-thermal emission associated with the ra- 
dio synchrotron halos observed in the halos 
of M82 and NGC 253 (c.f., Seaquist & Ode- 
gard 1991; Carilli et al. 1992). The GIS 
spectra are particularly vulnerable to these 
effects given the poorer spatial resolution and 
higher sensitivity above 2 keV of the GIS rel- 
ative to the SIS. However, we find that the 
counting rate in our background regions are 
consistent with spectra extracted from blank 
sky fields. Substituting the blank sky back- 
ground spectra for the background spectra de- 
scribed above does not significantly affect any 
of the spectral results discussed below. The 
telescope response was calculated assuming 
most of the counts within 6' are consistent 
with a point source (i.e., a source extent of < 
1'). We went to some considerable effort in 
investigating the effect of neglecting the spa- 
tial complexity on the spectral analysis and 
these results are described in the appendix. 
We find that neglecting the spatial complex- 
ity does not have a siginificant impact on our 
results. The spectra were binned to a mini- 
mum of 10 counts per bin to allow use of the 
X^ statistic. 

3. Spectral Fitting 

The SIS spectra strongly rule out a single- 
component model from 0.4-10.0 keV, with 
both power-law and Raymond-Smith plasma 
fits resulting in y^ > 2.0 for both sources. 
The poorer spectral resolution and response 
(at low energies) of the GIS allow for a single- 
component model. Indeed, most of the con- 
tribution to x^ i^i the SIS fits results from 
narrow features below ~ 2.5 keV. This can 
be seen in Figure 1 which shows the ratio of 
data to best-fitting power-law model, with the 
location of expected Kq, line emission from 
the He- and H-like ions of O, Ne, Mg, Si, 



Table 1 
Two-component Model Fits to ASCA M82 and NGC 253 Spectra 

Galaxy/Detector Model* A^^"-^* kT (keV) N^^'''^ T/kT (keV) A/A^q xV^of 



N253/S0-1 


1 


1 n+o-'^ 

-■-•U-i.o 


o.sitn^ 


15li? 


9 Q9+0.51 




416.9/360 


N253/S0-3 


1 


i.ot?:^ 


o.sotn^ 


141^3 


9 m+o-4o 




915.0/870 


N253/S0-1 


2 


o.sin 


n o 1+0.06 
U-8i_o.i6 


0+13 
^-8 


7 9+9.8 
'•^-3.2 


o.o3l°:i 


423.7/359 


N253/S0-3 


2 


n Q+0.9 
U-8_o.8 


o.sitH! 


81^^ 


7 Q+5.6 
1 ■o_2j 


o.o2t°:gi 


922.3/869 


M82/S0-lt 


1 


9 4+0-7 


0.64tSi^ 


1911^ 


1 7O+0.35 




548.1/422 


M82/S0-3 


1 


9 4+0-7 


0.63t°ig 


iq+io 


1 7f^+0.22 




1379/1334 


M82/S0-1 


2 


9 q+0.7 
^•'-'-2.3 


0.63;°;i 


i4_i4 


lOlf 


o.i9lS:t^ 


540.7/421 


M82/S0-3* 


2 


9 q+0.7 
^•<J-0.8 


0.63l°;i 


1/1+10 
i4_i4 


11.4111 


o.i8l°:i 


1366/1333 



Note. — The errors give the 90% confidence intervals assuming that all model parameters (ex- 
cept for normalizations) are interesting. Soft-component abundances derived from the 4-instrument 
Raymond-Smith plus power law fits are given in text. 

* 1: Raymond-Smith plus power law, 2: Raymond-Smith plus Raymond-Smith 
'^ Hard-component abundance relative to solar 

■*■ Due to a relative minimum in y^ space, the true soft component absoption column lower confi- 
dence limit may be zero (see text). 

iV^ is in units of 10^^ cm~^, applied to both continuum components 
]\j-hard jg ^^ uuits of 10^^ cm~^, applied to hard continuum component only 



S, Ca and Fe and the L-shell emission of Fe 
marked. The prominence of He-hke Mg, Si 
and S hnes imphes a thermal component with 
kT < 2 keV (c.f., Mewe, Gronenschild, & Van 
Den Oord 1985). Accordingly we tried fitting 
two-component models to the spectra, with 
a Raymond-Smith model to account for the 
soft (^ 2 keV) flux and a Raymond-Smith or 
power-law component for the hard {^ 2 keV) 
flux. These models resulted in acceptable fits 
for both NGC 253 and M82, although the SIS- 
only M82 fits resulted in x^ ~ 1.3 which can 
be rejected at a level of >99% (note that this 
neglects any systematic error). 

The results are shown in Table 1, where 
fits to the SIS data alone are denoted 'SO- 
1' and fits to the data from all four instru- 
ments are denoted 'SO-3'. Table 1 shows 
that the SO-1 and SO-3 fits give consistent re- 
sults. The errors given in Table 1 (and in 
the text below) assume that all model pa- 
rameters are interesting except for normal- 
izations, or 12 (8) interesting parameters for 
the M82 (NGC 253) Raymond-Smith plus 
power-law fits and 13 (9) interesting parame- 
ters for the double Raymond-Smith fits (the 
M82 and NGC 253 fits have a different num- 
ber of model parameters due to a difference 
in the handling of the soft-component abun- 
dances, as discussed below). The double 
Raymond-Smith model gives a better fit to 
the M82 spectra while the Raymond-Smith 
plus power-law model gives a better fit to the 
NGC 253 spectra, although the differences in 
X^ between the double Raymond-Smith and 
Raymond-Smith plus power-law models are 
statistically insignificant. The hard and soft 
components appear to be absorbed well in ex- 
cess of the Galactic columns (1.6 x 10^° cm~^ 
and 4.3 x lO^o cm^^ for NGC 253 and M82, 
respectively; Stark et al. 1992), although 
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Fig. 1. — Ratio of the SIS ASCA data to best- 
fitting power-law model for M82 (top) and 
NGC 253 (bottom). The possible identifica- 
tion of prominent lines is shown. 
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Fig. 2. — The best-fitting Raymond-Smith plus power-law model to the M82 and NGC 253 data. 
The top panels show the data with the model plotted as a solid line and the soft and hard 
components plotted separately as dashed lines, demonstrating that the hard component begins to 
dominate the spectra of both M82 and NGC 253 above ~ 2 keV. The SIS and GIS data are shown 
separately for clarity, with the SO and S2 points marked. The bottom panels show the ratio of 
the data to the models. As is evident in these figures, the SIS have a broader bandpass, with a 
low-energy cutoff of ~ 0.4 keV, than the GIS, with a low-energy cutoff of ~ 0.8 keV, while the 
GIS have more response at higher energies. 



not at a high statistical significance. The 
soft-component temperatures are similar (~ 
0.6 and 0.8 keV, for M82 and NGC 253, re- 
spectively, insensitive to the choice of hard- 
component continuum) while the hard com- 
ponent appears to be somewhat harder in 
M82 (r ~ 1.8 or kT - 11 keV) than in NGC 
253 (r ~ 2.0 or kT ~ 7 keV). If the hard com- 
ponent is thermal, the abundance (inferred 
from a Raymond-Smith model) is constrained 
to be significantly sub-solar (< 0.2 and < 0.5 
in NGC 253 and M82, respectively), based 
primarily on the lack of a strong Fe-K line 
at 6.7-6.9 keV. In NGC 253 no (narrow) Fe- 
K line emission is detected, and we set a 90% 
confidence upper-limit (for two additional pa- 
rameters, the energy and normalization of the 
line) of ~ 180 eV for the EW of any line emis- 
sion in the four-instrument Raymond-Smith 
plus power-law fit. In M82, the addition of 
a narrow line at ~ 6.6 keV with an EW of 
~ 100 eV to the four-instrument Raymond- 
Smith plus power-law fit reduces x^ by ~ 13, 
significant at a level of ~ 3(T for two addi- 
tional parameters. The Raymond-Smith plus 
power-law fits to the M82 and NGC 253 data 
are shown in Figure 2. 

In M82 the soft-component abundances 
(for the SO-3 Raymond-Smith plus power-law 
fit) relative to solar were found to be: N < 2.1, 
O = OAStlil Ne = 0.151°:?°, Mg = 0.22l°:lt, 



0,04+0-02 



Si = 0.35irA, S = OMlltr, Fe - w.w._o.o3- 
The Ar, Ca, and Ni abundances were tied 
together (due to limited signal-to-noise) re- 
sulting in a mean value of 0.23]'^o;23. Note 
that the Fe abundance inferred from the soft- 
component model is determined solely based 
on Fe-L emission, not Fe-K, due to the domi- 
nance of the hard component above ~ 2 keV. 
We note that these abundances were derived 
assuming that the soft component is entirely 



thermal. The poorer signal-to-noise of the 
NGC 253 data do not warrant allowing the 
elemental abundances to vary independently. 
O, Mg, Si, and S, which show strong line emis- 
sion in Figure 1, were grouped together, and 
Ne and Fe-L, which are not resolved from each 
other, were grouped separately. The resulting 
abundances were N, Ar, Ca, Ni < 0.66; O, Mg, 
Si, S = 0.27+[i:5I and Ne, Fe = 0.05i[J:Ji As 
in M82, these results are suggestive of a de- 
pression of Fe abundance relative to O, Mg, 
Si and S. 

The somewhat high value of x^ resulting 
from the SIS-only fits to the M82 data (see 
Table 1) suggests that a third component may 
be required. Relative minima in x^ space 
are present in the M82 fits when the soft- 
component temperature and N^ are varied. 
This results in a highly assymmetric confi- 
dence interval for the temperature (see Table 
1). In the case of Nh, the required difference 
in x^ is reached before the onset of the new 
minima, giving the lower limit shown in Ta- 
ble 1, however the "true" lower limit on N^ 
may be zero. This behavior also suggests the 
presence of an additional component, with a 
lower temperature than that of the soft com- 
ponent. The addition of a third thermal com- 
ponent to the double Raymond-Smith model 
substantially improves the M82 fits, resulting 
in temperatures of 0.35 keV, 1.07 keV and 
12.4 keV (see Tsuru et al, 1997 for a dis- 
cussion of three-component Raymond-Smith 
and non-equilibrium plasma fits to the M82 
spectra). Adding an additional thermal com- 
ponent to the double Raymond-Smith NGC 
253 model reduced x^ by only 8.3 for 2 addi- 
tional degrees of freedom (kT and normaliza- 
tion, with Nh fixed at the Galactic value and 
the abundances tied to the T ~ 10^ K compo- 
nent abundances), significant only at the 2.4cr 



level. The temperature of the additional com- 
ponent was ~ 0.14 keV and the best-fitting 
model parameters of the other two compo- 
nents were similar to the two-component fit 
values in Table 1. 

In Figure 3 we plot the abundances rela- 
tive to Fe (x/Fe) derived from various fits to 
the M82 data. In addition to the M82 fits in 
Table 1 and the triple Raymond-Smith model 
discussed above, these fits include substitut- 
ing the Mewe-Kaastra-Liedahl plasma model 
('mekal' in XSPEC) for the Raymond-Smith 
plasma model in the Raymond-Smith plus 
power law model, a power-law distribution 
of temperatures model, a double Raymond- 
Smith plus power-law model, and including 
an additional (very steep) power-law in the 
triple Raymond-Smith model. While these 
fits result in abundances that are one to two 
orders of magnitude supersolar (indicating that 
the soft-component continuum is poorly con- 
strained). Figure 3 shows that the relative 
abundances are somewhat more robust, with 
most varying by less than a factor of ~ 5, 
although the relative abundances of Ar, Ca 
and Ni are poorly constrained. Note that 
these errors are still somewhat larger than 
the statistical errors discussed above, show- 
ing that systematic errors dominate. These 
systematic errors are most likley dominated 
by differences in the slope of the continuum 
across the line-emitting energies in each fit. 
A comparison of the Raymond-Smith and 
the Mewe-Kaastra-Liedahl models in XSPEC 
for a 0.6 keV, solar- abundance plasma shows 
that, with the exeption of an Oxygen line 
at ~ 0.8 keV present in the Raymond-Smith 
model but not in the mekal model, most of the 
differences between the two models are on the 
order of 30% or less, implying that uncertain- 
ties in atomic physics are likely to be much 
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Fig. 3. — The abundances of various elements 
relative to Fe. The multiple points plotted 
for each element show the range of relative 
abundances derived from various models (see 
text). Note that the ranges shown are dom- 
inated by systematic errors which are much 
larger than the statistical errors {not explic- 
itly shown), such as those for the fits given in 
Table 1. Also shown are the yields predicted 
for 2OM0 and 25M0 supernovae, taken from 
Thielemann, Nomoto, & Hashimoto (1996), 
normalized to Fe. 



less significant than the statistical errors in 
these data. Also shown in Figure 3 are the 
x/Fe ratios expected for 20 and 25 Mq super- 
nova showing that even the direct emission of 
Type-II supernovae should not show such a 
large (relative) underabundance of Fe, partic- 
ularly relative to Si and S ( Thielemann et al. 
1996). 

In Table 2 we give the observed fluxes and 
absorption-corrected luminosities (given for 
each component separately) for NGC 253 and 
M82 in the 0.5-2.0 keV and 2-10 keV bands, 
with the SIS values corrected for non-uniform 
spatial exposure. The fluxes inferred from the 
individual detectors are within 3-4% of the 
mean values, although systematic differences 
of 10-20% between the SIS and GIS values are 
often found in other sources. The uncertainty 
in Nj:/ is signiflcant in computing the unab- 
sorbed luminosities, with the 90% confldence 
interval in N^^ in NGC 253 corresponding 
to an uncertainty of ~ 64% in intrinsic 0.5- 
2.0 keV luminosity. The 0.5-2.0 keV PSPC 
fluxes for M82 and NGC 253 from 6' regions 
are ~ 30% and 10% higher than the ASCA 
fluxes, respectively, and lie within the scat- 
ter of the correlation of the ASCA and PSPC 
fluxes for 13 low-activity galaxies (including 
M82 and NGC 253) shown in Serlemitsos, 
Ptak & Yaqoob (1996). Our 2-10 keV fluxes 
are ~ 1/2 the Ginga values (Ohashi et al. 
1990, Tsuru et al. 1990) although any inter- 
pretation of this result must be viewed with 
caution given the large (1° by 2°) Cinga FOV 
(however, the fact that the M82 observation 
was performed in a scanning mode reduces 
the uncertainty in flux). 

4. Discussion 



4.1. Comparison with Superwind Mod- 
els 

While two and three-component models 
describe the data well, the true situation, par- 
ticularly for the soft component below ~ 2 
keV, is likely to be a distribution of temper- 
atures. For example, the ROSAT PSPC data 
from both M82 and NGC 253 suggest tem- 
peratures of ~ 10^ K in their disk and nu- 
clear regions, with the temperature decreas- 
ing with radius out in the halos (i.e., beyond 
~ 1 kpc), as discussed in Strickland, Ponman, 
& Stevens (1996) and Dahlem, Heckman & 
Weaver (1996). In NGC 253, the soft (< 2 
keV) component is apparently dominated by 
a disk and nuclear thermal component with 
kT ~ 0.8 keV. In M82, on the other hand, 
the soft component appears to have approx- 
imately equal contributions from ~ 1.1 keV 
gas in the nuclear and disk regions and ~ 
0.4 keV gas in the halo, which is flt reason- 
ably well by a 0.8 keV thermal model. The 
luminosities, pressures, masses and energet- 
ics implied by the best-fltting soft- component 
model parameters are consistent with recent 
models of starburst-driven winds (Tomisaka 
& Bregman 1993; Suchkov et al. 1996), al- 
though the details of this kind of analysis are 
sensitive to assumptions concerning the vol- 
ume, fllling factor, and possibly the outflow 
velocity of the gas associated with each com- 
ponent. For example, the outflow velocity of 
the X-ray emitting gas in M82 may be as 
high as 1000-3000 km s"^ (Seaquist & Ode- 
gard 1991) although the outflow velocity of 
molecular gas in M82 is estimated to be less 
than 500 km s'^ (Nakai et al. 1987). To 
estimate the contribution of SN to the ob- 
served thermal flux we assume that the ISM 
heating is dominated by type-II SN, and that 
each supernova ejects ~ IOM0 with 10^^ ergs 
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Table 2: ASCA Fluxes and Luminosities for M82 and NGC 253 



Galaxy 


Energy Range 


Observed Flux 
(10~^^ ergs s"-"^ cm~^) 


Absorption-Corrected Luminosities 

(10^° ergs s~i) 


Two-comp. Fit^ 


Three-comp. Fit* 


Soft 


Hard 


Soft 


Medium 


Hard 


M82 

M82 
NGC 253 
NGC 253 


0.5-2.0 keV 
2.0-10.0 keV 
0.5-2.0 keV 
2.0-10.0 keV 


9.4 
20 

2.8 
4.7 


2.7 

0.2 

0.23 

0.03 


2.1 
3.5 
0.33 
0.37 


3.2 
0.0 
0.02 
0.00 


2.5 
0.6 
0.26 
0.03 


1.4 

3.4 

0.19 

0.35 



All values are averages from SO-3 fits [M82 (NGC 253) fiuxes from individual detectors are within 

3% (4%) of the mean]. 

Observed fluxes are inferred from the best-fltting two- component model (three-component models 

give similar results). 

^ Power law (hard) plus Raymond-Smith (soft) model 

* Triple Raymond-Smith model 



of kinetic energy (c.f., Suchkov et al. 1996). 
With these assumptions, independent of the 
assumed volume V and filling factor, 



M: 



SN 



[\nVkT + \mHnVvl^^^)M: 



10 



M, 



total 



likT 



lO^^rriHnV 



where Msn is the mass of gas ejected by SN 
(i.e., individual SN and the superwind), Mtotai 
is the total X-ray emitting mass, Mio is the 
mass ejected by a typical SN in units of lOM© 
and Viooo is the outflow velocity in units of 
1000 km s"^ We estimate (with Mio = 1 
and t;iooo = 1) that 4^ ~ 12 - 14% for the 



Mtotai 

1.1 keV and 0.8 keV gas (for M82 and NGC 
253, respectively) of the soft-component emis- 
sion. Mio and f looo are likely to be somewhat 
smaller than ~ 3, so this implies that a signifl- 
cant amount of the emission is due to ambient 
ISM, possibly "swept-up" by the superwind. 



as suggested by Suchkov et al. (1996). The 
low absolute abundances implied by some of 
the models of the soft component may there- 
fore be indicative of low abundances in the 
(ambient and swept-up) disk ISM. This cal- 
culation is not appropriate for the ~ 0.4 keV 
halo gas in M82 since energy dissipation is 
likely to signiflcant out in the halo. The un- 
usually high values of x/Fe shown in Figure 
3 indicate that the Fe abundance is "sup- 
pressed", particularly relative to S and Si. 
The high IR luminosities of M82 and NGC 
253 suggest a signiflcant dust distribution in 
the disk and halo ISM of M82 and NGC 253 
(Telesco 1988). Since Fe is depleted in warm 
ISM clouds to a larger extent than S and Si 
(c.f., Cowie & Songaila 1986), dust depletion 
may account for some of the Fe deflciency, re- 
sulting in Si/Fe and S/Fe ratios ~ 17% and 
130% (respectively) higher than in the unde- 
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pleted case 0. 

Repeating the "swept-up mass" calculation 
described above for the hard component, we 
find that, if thermal, ~ 25-75% (NGC 253) 
and 40-100% (M82) of the hard component 
originates in superwind emission and not the 
ambient ISM. The high absorption {Nh ~ 
10^^ cm~^) derived from fits to the hard com- 
ponents (see Table 1) also suggests that the 
hard component originates in compact and/or 
nuclear regions (note that radio observations 
of the nuclear regions of both M82 and NGC 
253 suggest the presence of molecular gas with 
similar column densities; Nakai et al. 1986; 
Paglione, Tosaki, & Jackson 1996). How- 
ever, the inconsistency of the low abundances 
inferred from the hard- component fits with 
the high metallicity expected in supernova 
ejecta and solar abundances expected in the 
nuclear HII regions of M82 and NGC 253 
(c.f., Carral et al. 1993; Lord et al. 1996) 
suggest that the hard components may not 
be dominated by (nuclear) superwind emis- 
sion. Dust depletion may again be impor- 
tant although dust destruction is likely to be 
significant in nuclear starburst regions over 
the ~ 10^~^ year lifetime of the starburst ( 
Calzetti, Kinney & Storchi-Bergmann 1996). 
It is therefore likely that gas not in ther- 
mal equilibrium or non-thermal emission con- 
tribute significantly to the hard components. 
The former case would correspond to uq < 
0.5 - 1.5 X 10~^(i;iooo/'"fcpc) cm"^, where rupc 
is the hard-component emitting region size in 
kpc, f 1000 is the velocity of the ejecta (as de- 
fined above), and uq is the pre-shock ambi- 
ent density (c.f., Tomisaka & Bregman 1993; 



^The dust depletion given in table 3 of Cowie & 
Songalia was scaled to the '^ 50% depletion of Si es- 
timated by Lord et al. (1996) for the nuclear HII 
regions in M82 in this calculation. 



Hamilton, Sarazin, & Chevalier 1983). Even 
if "^1000 is as high as 3, such a low value of 
riQ is unlikely in the nuclear regions of M82 
and NGC 253 (Carral et al. 1993; Lord et 
al. 1996), suggesting that non-thermal equi- 
librium ionization is not important unless the 
direct emission of SN is significant. Using the 
radio surface brightness - radius relationship 
for SN in Huang et al. (1993) and the X- 
ray luminosities of SN1986J and Gas A as a 
guide, we expect the X-ray luminosity of SN 
in M82 to be ~ lo^^-ss ^^^^ g-i^ -j-qq j^^ ^q 

account for more than a few percent of the 
hard component, while in NGC 253 some of 
the bright extranuclear sources may be super- 
novae contributing to the hard component. 
Future work will examine the possible ther- 
mal contribution to the hard component in 
more detail, but it appears likely that the 
hard components in M82 and NGC 253 are 
dominated by non-thermal emission. 

4.2. The Nature of the Non-thermal 
Emission 

The two most likely candidates for non- 
thermal emission are accretion-driven point 
sources (e.g.. X-ray binaries and AGN) and 
inverse- Compton (IC) scattering of (nuclear) 
IR fiux by relativistic electrons. The IC 
component may be significant (estimated by 
Schaaf et al. 1991 to be Li.4_8.9 kev ~ 1-5 x 
Xq40 gj,gg g-i Qj^ ^ bQ% of the hard compo- 
nent in M82), although the likely high ratio 
of proton to electron energies (i.e., ~ 100; 
Seaquist & Odegard 1991) would significantly 
reduce the contribution of this component. 
On the other hand, the point sources in M82 
and NGC 253 may account for much of the 
hard component. In M82, a nuclear source 
has varied from ~ 0.04-0.06 counts s~^ to ~ 
0.02 counts s^^ between the 1991 ( Collura et 
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Fig. 4.— The ROSAT HRI background- 
subtracted light curve for the nuclear source 
in M82 (c.f., Figure 2 of Collura et al. 1994), 
including recent AO-5 data. The source 
counts were extracted from a 10" region 
centered on the source and the background 
was determined from a 15" -20" annulus cen- 
tered on the source (here the background 
is completely dominated by nuclear emission 
in M82). Also shown is the corresponding 
absorption-corrected 2-10 keV luminosity as- 
suming a spectrum similar to the hard com- 
ponent inferred from the ASCA data (i.e., 
r ~ 1.8 and A^^ ~ 10^^ cm"^). Evidently, 
this source could account for a considerable 
fraction of the intrinsic 2-10 keV luminosity 
of M82 {Lx ~ 3.4 x 10^° ergs s^^) observed 
by ASCA (on day 755 in this figure). 



al. 1994) and 1995 ROSAT HRI observations, 
corresponding to an intrinsic 2-10 keV lumi- 
nosity of l-4x 10^° ergs s~^ or ~ 30%-100% of 
the hard- component luminosity assuming this 
source has a spectrum similar to the power- 
law model fit to the hard component (see Ta- 
ble 1 and Figure 4). In NGC 253, four strong 
extranuclear point sources (~ 2-4' from the 
nucleus) are responsible for ~ 20% of the 0.5- 
2.0 keV flux in the ROSAT PSPC observa- 
tion. As shown in Figure 5 (plate 1), the point 
sources in NGC 253 have varied between the 
Einstein and ROSAT observations and are 
clearly prominent above 2 keV in the ASCA 
SO image. We extracted composite SIS spec- 
tra for the two brightest extranuclear sources 
(sources Rl and R5 in flgure 5), using source 
region sizes of 45" to minimize the "contam- 
ination" from diffuse gas. Fitting these spec- 
tra with the Raymond-Smith plus power-law 
model (see Table 1) resulted in fit parameters 
consistent with those listed in Table 1 for the 



2.1+1 X 10 



^1.9 



22 



cm 



entire galaxy (e.g., Nh 

and F ~ 2.3l'^|;o). These fits suggest that 
sources Rl and 5 account for ~ 25% of the 
observed 2-10 keV fiux from NGC 253. While 
it is possible that some of the point sources in 
NGC 253 are associated with SN (see above), 
some are likely to be X-ray binaries given the 
lack of optical counterparts (i.e., Fabbiano & 
Trinchieri 1984). The Eddington luminosities 
of these sources suggest masses of at least 3-20 
Mq implying that they are blackhole candi- 
dates, which also have spectra consistent with 
the hard components observed in M82 and 
NGC 253 (Tanaka 1989). The spectra of M82 
and NGC 253 are also very similar to the spec- 
tra of low-luminosity AGN and LINERs (Ser- 
lemitsos, Ptak and Yaqoob 1996) which gen- 
erally can be described by a Raymond-Smith 
plus power-law model. The power-law slopes 
obtained from our fits (F ~ 1.8 — 2.2) may 
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suggest a common physical origin, such as an 
accretion-driven component. It is hkely that 
some of the point sources in M82 and NGC 
253 have masses ^ IOMq since it would be 
unphysical for all of these sources to be emit- 
ting at their Eddington luminosities most of 
the time, if they are powered by accretion. 
The discovery of (long-term) X-ray variabil- 
ity in "starburst" galaxy NGC 3628 demon- 
strates that the X-ray flux from NGC 3628 
is dominated by a point source, most likely a 
massive blackhole candidate or a micro- AGN 
( Dahlem, Heckman, Fabbiano 1995; Yaqoob 
et al. 1995). These considerations suggest 
the possibility that starburst galaxies such 
as M82 and NGC 253 may harbor blackhole 
candidates, with masses significantly higher 
than the masses of blackhole candidates in our 
galaxy. 
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Appendix: Computation of the ASCA 
X-Ray Telescope (XRT) Response for 
Extended Sources 

As stated in the text, the X-ray emission 
in both M82 and NGC 253 is extended while 
the telescope response was computed assum- 



ing it is point-like. In this appendix we show 
that properly taking the extended nature of 
the sources into account does not significantly 
impact the results of this paper. In X-ray 
spectral fitting, the XRT response is taken 
into account by the creation of a table of ef- 
fective area as a function of energy, which is 
then used to convert the observed counting 
rate in each spectral bin to a flux. This table, 
hereafter referred to as an ARE (ancillary re- 
sponse file), also applies a correction for the 
counts scattered out of the source region (6' 
circular regions were used for both the SIS 
and CIS in this paper). 

The ASCA XRT Response for a Point 
Source 

Eirst, we describe the "standard" proce- 
dure used for the creation of an ARE ap- 
propriate for a point source. The telescope 
point-spread function (PSE) is summed over 
the source region to determine percentage of 
source photons contained in the source re- 
gion. The PSE sum is then multiplied by the 
telescope effective area. The PSE and effec- 
tive area are functions of both energy and the 
off-axis angle of the source. Both the PSE 
and effective area are interpolated from values 
contained in calibration look-up tables, com- 
puted using a Monte-Carlo procedure based 
on a raytracing program. Although the XRT 
PSE has a cross-like shape, the calibration ta- 
bles contain azimuthally-averaged values, an 
approach that is usually a good approxima- 
tion as long as most of the source counts lie 
in a symmetric (and preferably circular) re- 
gion. The mirrors contain a heat shield, made 
up of a thin plastic film, to protect against 
temperature gradients, and the computed ef- 
fective area is reduced to take photon losses 
through the shield into account. The gas cell 
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Fig. 5.— Einstein HRI (0.2-4.0 keV, 7/5-9/79; top left), ROSAT HRI (0.1-2.2 keV; 6/7/92; 
top right), ROSAT VS^VG (0.5-2.0 keV; 12/6/91-5/6/92; bottom left), and ASCA SO (2-10 keV; 
5/25/93; bottom right) images of NGC 253. Each image is 10' by 10' and has been adaptively 
smoothed to produce a SNR of ;^ 5 at each pixel (without background subtraction). The sources 
discussed in Fabbiano & Trinchieri (1984) are marked (sources 3 and 6 lie outside of the FOV of 
these images) along with the X-ray nucleus ('N') and the "new" i?05'j4r source ('Rl'). 
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of the GIS contains a Beryllium window held 
in place by a supporting grid, and transmis- 
sion through the grid and Be must also be 
taken into account. Note that the thicknesses 
of both the grid and Be vary across the GIS 
detectors. Also, unlike the SIS, the GIS in- 
troduce a broadening of the PSF. The in- 
strumental PSF of the GIS is approximately 
a Gaussian whose FWHM ~ O.S^^^J^ ar- 
cminutes, with the addition of a "tail" at en- 
ergies ^ 4 keV. Accordingly, the XRT PSF is 
modified by the grid and Be transmission and 
then convolved with the GIS PSF, which is 
modeled by the sum of two-dimensional Gaus- 
sians. Hereafter, we refer to ARFs computed 
with this approach as standard ARFs. 

Computation of the ASCA Response to 
an Extended Source 

There are two ways in which standard 
ARFs may be inappropriate for an extended 
source. First, more source counts would be 
scattered out of the source region than is the 
case for a point source. Not taking this into 
account would result in an underestimate of 
the source flux. Second, by definition, an ex- 
tended source subtends a distribution of off- 
axis angles, so determining the XRT effec- 
tive area and PSF for a single off-axis angle 
may be incorrect. Because both the effective 
area and PSF have an energy dependence, the 
spectral results may be affected as well as the 
flux estimates. Because the X-ray emission 
from M82 and NGC 253 is relatively compact 
and weak, it is likely that statistical errors will 
dominate that the effects described above, 
with the possible exception that the fluxes 
may be underestimated. To check this in de- 
tail, we adopted a Monte-Carlo approach. As 
is the case in the point-source procedure, we 
divided the bandpasses of the SIS and GIS 



into the same number of energy bins as is 
used in the spectral fltting, namely 1180 in 
the case of the SIS and 201 in the case of 
the GIS. For each bin at least 10000 photons 
were raytraced through the mirror (the ray- 
tracing code was developed by Peter Serlemit- 
sos during the design of the ASCA mirrors, 
and, when azimuthally-averaged, is in good 
agreement with the look-up tables described 
above and the radial proflle of an ASCA PV 
phase observation of 3C 273). The position 
(on the sky) of each photon was chosen at 
random from ROSAT PSPC images of M82 
and NGC 253, weighted by the intensity at 
each pixel. This approach results in photons 
incident on the XRT with a distribution in 
off-axis angles consistent with those expected 
from the PSPC images. A 12' by 12' im- 
age was used for M82 and a 10' by 10' im- 
age was used for NGC 253. The bandpass 
of each image was limited to ~ 0.5 — 2.0 keV, 
and the background was not subtracted, in or- 
der to conservatively overestimate the source 
extent (note that the source extent is also 
overestimated since the ROSAT PSPC PSF 
is not removed with this procedure, although 
the ROSAT FSF FWHM is ~ 10", somewhat 
smaller than the ASCA XRT PSF FWHM of 
~ 50"). 

In the case of the GIS, the instrumental 
PSF was introduced by adding "noise" (i.e., 
Gaussian random deviates) to the raytraced 
photon positions. The standard deviations of 
the deviates were computed with the same 
algorithm used in the standard ARF proce- 
dure. Also in the case of the GIS, photons 
were rejected to account for the grid and Be 
transmission (i.e., a photon was rejected if a 
uniform deviate exceeded the product of the 
Be and grid transmission at the photon's posi- 
tion). The percentage of photons "detected" 
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in the source region, reduced by the thermal 
shield transmission (as is done in the point 
source case), was then multiplied by the geo- 
metrical area of the mirror to derive an effec- 
tive area. Hereafter, we refer to ARFs com- 
puted with this procedure as "raytraced" (al- 
though, as stated above, the "standard" ARF 
procedure also indirectly involves raytracing, 
in the computation of the look-up tables). 

Initially, ARFs were raytraced assuming 
that M82 and NGC 253 are point sources for 
comparison with the standard ARFs used in 
the fits given in the text. The ratio of the ray- 
traced ARFs to the standards ARFs, shown 
in the top panel of figure 6 for the case of 
the M82 SO and S2 ARFs, gives the system- 
atic differences between the two procedures 
described above. The largest source of dis- 
crepancy between the two methods is due to 
differences in the raytracing code used here 
and the raytracing used to produce the look- 
up tables, particularly in the choice of opti- 
cal constants for gold (the reflecting material 
used in the XRT). Other sources of discrep- 
ancy are inherent to the procedures employed, 
i.e., the XRT PSF is not azimuthally-averaged 
in the raytraced case. The top panel of figure 
6 shows that, with the exception of a sharp 
peak around the Au M edge at 2.2 keV, the 
two procedures are in agreement below 4 keV 
to within 10%. Above 4 keV there is a broad 
"hump" in the ratio plot which peaks at a 
level of ~ 20%. These are issues which are 
currently being resolved by the instrument 
teams and will be reported in a future issue 
of ASCANews. For comparison, the statisti- 
cal errors in the ASCA data below 4 keV are 
typically on the order of 6% in the case of 
M82 and 10% in the case of NGC 253. As is 
evident in Figure 2, the statistical errors be- 
come somewhat larger above 4 keV since both 
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Fig. 6. — (top) The ratio of ARFs for a point 
source computed using the raytracing algo- 
rithm described in the text and the "stan- 
dard" procedure using look-up tables for the 
XRT effective area and PSF. The ratio for 
ARFs appropriate for the SO (solid line) and 
G2 (dashed-line) M82 data are shown. The 
largest source of discrepancy between the ray- 
traced and "standard" ARFs arises from dif- 
ferences in the raytracing code used in this 
analysis and the raytracing used to produce 
the look-up tables, {bottom) The ratio of 
the ARFs computed assuming that M82 is 
extended and assuming that M82 is a point 
source in the case of SO (solid-line) and S2 
(dashed line). Here, both the point-source 
and extended-source ARFs were raytraced to 
remove the discrepancies between the ray- 
traced and "standard" procedures (shown in 
the top panel). 
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the XRT and instrumental responses and the 
source fluxes drop off rapidly above 4 keV, 
particularly in the case of NGC 253 which has 
a steeper spectrum than M82. 

Implications for M82 and NGC 253 

The bottom panel of figure 6 shows that 
ratio of ARFs produced with the raytracing 
procedure for an extended source (with pho- 
tons selected from the PSPC image, as de- 
scribed above) and a point-source (raytraced 
ARFs were computed in both cases to remove 
the discrepancy between the raytraced and 
standard ARFs) in the case of the M82 SISO 
(solid line) and GIS2 (dashed line) ARFs. 
Note that the extended nature of the source 
only impacts the overall normalization be- 
low ~ 4 keV, and results in a gradual drop 
off in net effective area above 4 keV (a re- 
sult of the larger energy dependence of the 
XRT PSF and effective area above 4 keV). 
Note that this drop in effective area par- 
tially offsets the "hump" above 4 keV noted 
above. Indeed, in the case of NGC 253, re- 
peating the four-instrument Raymond-Smith 
plus power-law fit with the extended source 
ARFs results in best-fitting parameters that 
are essentially identical to those discussed in 
the text where the standard ARFs were used 
(i.e., the statistical errors dominate any sys- 
tematic errors). When the double Raymond- 
Smith fit to the M82 data is repeated with 
the extended-source ARFs, the soft compo- 
nent parameters do not change significantly, 
although the inferred 0.5-2.0 keV fiux is in- 
creased by ~ 40%. Because our procedure 
dehberately overestimates the source extent, 
this is a conservative upper limit to the 0.5- 
2.0 keV fiux from M82. On the other hand, 
the hard component spectral parameters are 
modified slightly when extended-source ARFs 



are used, with the best-fitting N^ decreasing 
by ~ 5 X 10^^ cm^^, the temperature increas- 
ing by ~ 3 keV and the abundance increasing 
from ~ 0.2 solar to ~ 0.3 solar. Nevertheless, 
these parameters are still well within the sta- 
tistical errors quoted in the text. As with the 
soft component, the extended-source fits sug- 
gest that the 2-10 keV fiux, dominated by the 
hard component, may be underestimated by 
~ 27%. However, Tsuru et al. (1997) report 
that the fiux associated with the hard compo- 
nent in M82 is not resolved by ASCA. Recom- 
puting the extended-source ARFs by selecting 
the soft-component photons from the PSPC 
image and assuming that the hard compo- 
nent is point- like eliminates the changes in the 
hard- component fit parameters, and results in 
a 2-10 keV fiux ~ 7% lower i\ia,n that inferred 
from the standard ARF fits (this is not sur- 
prising since, as shown in the top panel of Fig- 
ure 6, the raytraced point-source ARF gives a 
larger effective area at high energies than the 
standard ARFs). Thus, the only parameter 
significantly affected by properly computing 
the response of ASCA to an extended source 
is the 0.5-2.0 keV flux inferred for M82, which 
does not impact any of the conclusions of this 
paper. 
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